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Abstract 

We study the algebraic structure of an integrable Hubbard-Shastry type lattice model associated with the centrally 
extended su(2|2) superalgebra. This superalgebra underlies Beisert’s AdS/CFT worldsheet R-matrix and Shastry’s 
R-matrix. The considered model specializes to the one-dimensional Hubbard model in a certain limit. We demonstrate 
that Yangian symmetries of the R-matrix specialize to the Yangian symmetry of the Hubbard model found by Korepin 
and Uglov. Moreover, we show that the Hubbard model Hamiltonian has an algebraic interpretation as the so-called 
secret symmetry. We also discuss Yangian symmetries of the A and B models introduced by Frolov and Quinn. 
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1. Introduction 

Exactly solvable models of strongly correlated electrons 
are of great importance in theoretical condensed matter 
physics. For instance, they play a prominent role in un¬ 
derstanding high-Tc superconductivity. The key example 
of such a model is the one-dimensional Hubbard model 
introduced in [1] . It describes the dynamics of interacting 
electrons in a one-dimensional lattice that models the con¬ 
duction band of a solid. Each site in the lattice can have 
four different states. It can be vacant, occupied by a spin 
up or down electron or occupied by an electron pair. 

The Hubbard Hamiltonian Jf = -|- 

for the inhnite lattice is written in terms of the usual 
creation and annihilation operators c^a, with a = f, 
of electrons with spin up and down. Here 

^ii+l •= 'y ^ Vj := (Cj^— 2)’ 

where tiic := cj^Cia is the number operator, i G N enumer¬ 
ates lattice sites, ^ := and h is the coupling constant. 

The integrability of the Hubbard model is known since 
the works of Lieb and Wu [2]. However, the R-matrix was 
only found much later by Shastry [3]. Given an integrable 
model, the R-matrix can usually be found from the under¬ 
lying symmetries of the model. Such symmetries of the 
Hubbard model were unknown until very recently. 
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The Hubbard model was long known to have an exact 
so(4) = (su(2) 0 su( 2 ))/Z 2 symmetry [HIS] that can be 
extended to a certain Yangian [5]. Nevertheless, this was 
not sufficient to determine Shastry’s R-matrix and it was 
suspected that there are more symmetries underlying the 
Hubbard model. 

The answer to this question of the full symmetry al¬ 
gebra of the Hubbard model came from an unexpected 
area - superstring theory. The R-matrix of the Hubbard 
model emerges in the prime example of the AdS/CFT cor¬ 
respondence as the worldsheet R-matrix laiHiE]- It is 
the intertwining matrix for fundamental representations 
of the centrally extended su(2|2) superalgebra. We will 
denote this superalgebra by g. Interestingly, this superalge¬ 
bra has a non-standard Yangian extension, which is also a 
symmetry of the R-matrix [10) . An exceptional feature of 
this Yangian is that it has an additional generator, which 
has no counterpart in g, called the secret symmetry HU. 
The complete Yangian symmetry is an infinite-dimensional 
superalgebra of a novel type [12) . 

In this letter we consider an integrable one-dimensional 
lattice model by identifying each site of the lattice with a 
fundamental module of g. This gives an integrable Hubbard- 
Shastry type model which specializes to the Hubbard model 
in a certain limit. Such generalized models were first con¬ 
sidered in m and El- We call this model the general 
model. Likewise, there are three other parity invariant 
lattice models that follow from the general model, - the A 
and B models of Frolov-Quinn m and the Essler-Korepin- 
Shoutens (EKS) model [TO] . 

We construct an oscillator representation of g, its Yan- 
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gian extension and the R-niatrix, and discuss the integrable 
structure of the general model and its specializations. In 
particular we shed more light on the symmetries of the 
Hubbard model and obtain the spin Yangian symmetry 
of the Hubbard model [ 6 ] as a certain specialization of 
the Yangian of g. Moreover, we show that the Hubbard 
Hamiltonian is a specialization of the secret symmetry. We 
also obtain Yangian symmetries of the A and B models. 

This letter is organized as follows. We first introduce 
the necessary notation, the setup of the lattice, and the 
algebra g. We then construct an oscillator representation 
of g and its Yangian extension, which allows us to compute 
the R-, Lax- and transfer matrices of the general model. 
In the remaining sections we discuss links with the A, B 
and EKS models, the Hubbard model limit, and prove the 
aforementioned results. 


where 1 = 4 = 0 and 2 = 3 = 1 is the grading. A matrix 
representation of the oscillator algebra is given by the map 
TTi : OsCi —End(I/i) so that 

c*; ^ - E 34 ), Qt ^ + E 24 ), 

C4 ^ - E43), ^ e-^^^\E 3 i + E42), 

^4 I— E 22 + E 44 , 04 I—)■ E 33 + E 44 . (2.2) 

The additional phases (j)a in the expressions above are 
required to make the Hamiltonian and the R-matrix of the 
model manifestly invariant (i.e. without additional twists 
or similarity transformations) under the standard spin and 
charge symmetries. This will be explained in detail in 
Section 4. For subsequent reference we also set (-)'*' := 1 
and (-)^ := - 1 , (-)tl := - 1 . 


2. Oscillators and vector spaces 

Oscillator algebra. We consider a one-dimensional lattice 
with L e N sites that can be occupied by spin-^ particles. 
Each lattice site is identified with a four-dimensional Z 2 - 
graded vector space Vi = spanned by vectors 

Vi = spanc{ |0),, 14)*, |t)i, |t4)i }, 

where 1 < i < L is the index of the site in the lattice. The 
entire lattice is the L-fold tensor product V := 0^^ Vi. 
The vector |0)i is the vacuum (i.e. a hole at the site), |t)i 
and | 4 )i are spin up and spin down particles, and |t 4 )i is a 
multi-particle state occupied by a pair of spin up and spin 
down particles. The grading is 1 for vectors |t)* and |4)i, 
and 0 otherwise. To describe such a lattice, we introduce 
the usual free-fermion oscillator algebra generated by 
and Cia with 1 < i,j < L and a,/3 = t)4i that satisfy the 
standard (anti)commutation relations 

{CionCjp\ = 0 , = 0 , = \5ij5afi. 

We will denote the universal enveloping algebra of a single 
copy the oscillator algebra by Osc^. Likewise, 0 sC 4 (fc,..) 
will denote the universal enveloping algebra oiij{k...) copies 
of the oscillator algebra. For ease of notation, we identify 
the unit 1 of the oscillator algebra with the unit 1 of the 
ground field C. Upon identification 

|a)* := e^^-^^cLlO)*, |t4). := ( 2 . 1 ) 

and requiring CiQ,|0)i = 0, the space Vi becomes a left Osc^- 
module; here i := and (j>a & V. are arbitrary phases. 

Matrix representation. Let Eab S End(Ui), with 1 < a, 6 < 
4, be the usual graded matrix units satisfying 

[Eab,E,4 = SbcEad - 

(Eab ® E,d)iE,, 0 Eki) = EabE,, 0 E,dEki, 


3. Superalgebra 

Superalgebra. We recall the definition of g due to m By 
[•, •] we will denote the graded commutator [a, b] = ab — 
{—l)°‘^ba for any elements a and b in the superalgebra. The 
bar “ : g —Z 2 denotes the degree of the element under 
the Z 2 -grading. We will call grade 0 elements bosonic. 
Likewise, grade 1 elements will be called fermionic. 

The centrally extended superalgebra su(2|2) is gener¬ 
ated by elements Ea, Ea, Ha, with a = 1 , 2 ,3, and central 
elements P and K subject to the following defining rela¬ 
tions: 

[Ha, Eb] = AabEb, [Hi, Ej] = —AabEb, 

[Ea,Eb] = SabDaHa, [E 4 , E 3 ] = 0, [^ 1 ,^ 3 ] = 0, 

[Ea,[Ea,E 2 ]]^ 0 , [Fa,[Fa,E 2 ]]= 0 , 

[[E 4 ,E 2 ], [E 3 , E 2 ]] = P, [[Fi, F 2 ], [F 3 , F 2 ]] = K 

for a,b = 1,2,3 and c = 1,3. Here D = diag(l, —1, —1) 
is the normalization matrix and A is the Cartan matrix 
given in the appendix. Dynkin nodes 1 and 3 are bosonic; 
Dynkin node 2 is fermionic. Consequently, the Z 2 -grading 
is 1 for E 2 , F 2 and is 0 otherwise. 


Local oscillator representation. Fix a complex number h 
(this will be the coupling constant of the model). To each 
site we associate a spectral parameter Ui, so that the lattice 
is parametrized by u = (ui,U 2 , and h. 

To obtain a representation of g we introduce a function 
x{u) and parameters defined by m 


x{ui) := - 


Ui 


-4 


:= x{ui ± 4^)* . (3.2) 


The parameters xf^ are Zhukovsky variables satisfying 

-el _ 1 , + 1 _ 1 „ 

xj — 1 : - X^ - - = n, xj A —v + iCj -I— 3 = 2ui. 

Xj X, Xj X: 


Then, for each site, we introduce the local weights 


_ T'* e _ “ 1 _ J _ 

“ /fc ’ “ /fc ^ ~ PR ~ 

Vn v«7i VnxJ 


- 1) 
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where vf = xf lx~ and = Vi{xf — x~). Local weights 
describe the fundamental representation of g at each siteQ 
A local oscillator representation of g is given by the 
map : g —>■ OsCi such that (for a = 1, 2,3) 


^ ^i,a: Fa ' ^ I ^ 

K ^ Xi, C ^ Cl, 


(3.3) 


where 

£..i = = e-*(^t+<^r)*C4Qt, 

£*.3 = J,.3 = e*(^^-^t)*ct^c4, 

£*.2 = e“*‘^^*((ai - &i)c4n4 + 6iC4), Ti = aA, (3.4) 
3^i,2 = e {{d-i ~ Ci) Cj^n 4 + c^^), X^ = Cidi, 

= n4+ n4 —1, = n4 — n4, 

dfi,2 = —Ci — Ci = \{aidi + biCi) 


are operators in Osc^ C OsCi.../,- The operators with labels 
1 and 3 form two copies of the su(2) algebra. At the A 
site, for a = f, 4-, we have 


£^,3|i)^ = It)*, T,,3|t)* = It)*, df,,3|a)i = (-)“|a),. 

Thus operators with label 1 realize the spin symmetry 
of the spin-i particles. Similarly, operators with label 3 
satisfy 

£*.i|0)* = Itt)*, T,,i|tt)* = |0)„ Jf,,i|a), = -(-)“|a), 

for a = 0,tt- This is the charge symmetry of the spin-^ 
particles. From now on we will refer to these operators 
using this terminology. The action of fermionic operators 
(for a = t, tt and /3 = 0, f) is 

£*,2|t)* — ^*|0)*, £*,2|tl')* — a^lt)*, dfj^lo)^ — 2^*^*!^)*, 

Ti,2|0)i = Ci|t)i, Ti,2|t)i = rfi|tt)i, dfi4|/3)i = -\biCi\P)i. 
These operators are supersymmetries of the model. 


Global oscillator representation. Let us extend the repre¬ 
sentation (3.41 to the entire lattice. To achieve this, we 
need the tensor product describing the lattice V = {S>^i ^ 
to be defined over the oscillator algebra OsCi...^. We re¬ 
quire 3 i,a \a)j =±\a)j 3i,a if * j for all 3 i,a in ( [0| ) and 
a S Here the sign is ’ if both 3i,a and |a)j 

are fermionic, and ‘-I-’ otherwise. Following the braided 
Hopf algebra structure of g introduced in [HI, we define 
the so-called braiding factors = ni<j<* set 

iy(o) ._ g^j^d (for 1 < a < 3) 


l<i<L 


£a= E T=E(^ 

l<i<L 

Ta= ^ X = 

l<i<L 

dfo = ^ ^i,a, 






e= ^ e,. ( 3 . 5 ) 


l<i<L 


l<i<L 


l<i<L 


^The fundamental representation of g is a four-dimensional atypical 
cyclic (non-highest weight) representation. See sections 2.3 and 2.4 
in [7| for more details. 


Notice that the braiding factors make the action 

of the fermionic operators £ 2 , T 2 and the central elements 
T, D<C non-local. We will call operators in (3.51 the global 
operators. 

The map p ; g —>■ OsCi...l obtained from the one in 
(3.3) by replacing the local operators in (3.4) with their 
global couterparts in (3.5) is a global representation of g. It 
is an oscillator algebra realization of the braided coproduct 
introduced in m- 

For any Cja,cJ„ e OsCi.,,l let 


(j : 



^ i 0«CL-22+1) t 

^ ^ ^L-z+l,a 


(3.6) 


be the parity reversal operator. Moreover, require 17 ( 0 ^) = 
for all local parameters (i.e. local weights, Zhukovsky 
variables, etc.). Then the ‘opposite’ map p' := a o p is a 
realization of the opposite braided coproduct. 


Yangian. The superalgebra g admits a non-conventional 
Yangian extension Y(g) [1^. To describe its structure 
we need to introduce local operators Qiap and corre¬ 
sponding to all fermionic root vectors of g: 


^*tt — £*, 2 , 

= [£*, 2 , £*,i], 
^*T4 = [£*,3, £*, 2 ], 
^44 = [£*,i, Qip], 


^44 - 

= [^M,^*. 2 ], 

Qin = [^^2,T*.3], 

Q !44 = [Q 1 , 2 ,^.i]- 


(3.7) 


Their global counterparts, and are defined equiv¬ 
alently. In addition, we define the following non-local 
quadratic operators 


2,4 — 'y ( Da{£-i^a3^j,a + 3^i,a£j,a + ■, 

a=l,3 


a,/3=ti 


yO-l) ^ ^ 

y{i-l) =*= 4-1) 2*0/3 


For 3 € {£, T, df} we then introduce spin and charge Yan¬ 
gian operators (for a = 1, 3) 

3a = [3i,aAij — (3.8) 

i<j<L 


The Yangian operators corresponding to the remaining 
generators of g can be obtained using Lie algebra relations. 
Finally, we introduce an additional Yangian operator 


■B = 


y —^ 


djCj 




I E (3-9) 

1 < 2 <J'<L 


The operator 23 is called the secret symmetry m- It has 
no counterpart in the algebra g and extends the Yangian 
of g to a deformed Yangian of g[(2|2) [12]. 
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PPHT. We introduce a partial particle-hole transforma¬ 
tion P^, which acts as the identity map on the spin up 
oscillators, but as involution on the spin down oscillators, 
and inverts the lattice parameters as follows: 

Pi ■ ^ (3.10) 

This transformation interchanges the spin and charge oper¬ 
ators, £ip O —£i, 3 , Tip O —Tip, TCip o —IKip; inverts 
the elements Ti ——Ti, 3Ci —>■ —3Ci, Ci —>■ —Ci; and acts 
on the fermionic operators as £ip i—>■ ^£i, 2 , Tip t ^Tip 
due to the presence of Vh factors in the local weights. 


4. Integrable structure 


R-matrix. The R-matrix of an integrable system is de¬ 
termined by its symmetry properties. For a subset S C 
{1, ...,T} introduce a projector es : OsCi..,l — t Oscg by 

es(cia) = es(cL) = ^iescL, es{i^i) = (^'i)'^'^®, 


where Si^s = 1 if i G T, and Si^s = 0 otherwise. Then, 
for each of the global operators 0 defined in (3.5) and also 
for the Yangian operators in (3.8) and (3.9), the oscillator 
R-matrix is the two-site operator 5lij G OsCij satisfying 
the intertwining equation 


Tij eij{d) = {eij o a){3) (4.1) 

This equation implies that the oscillator R-matrix is a 
homogeneous element of degree 0 in OsCij, and determines 
it uniquely up to an overall scalar factor (which we set to 
identity for simplicity). 


where It is an instructive exercise to 

check this oscillator R-matrix is unitary, JiijJlji = 1 , and 
satishes the Yang-Baxter equation 

= 3ljk3lik3lij. (4.4) 

The transformation P^ is not a manifest symmetry of Jlij. 


Matrix representation. The graded matrix representation 
R{ui,Uj) of the oscillator R-matrix is obtained by setting 

R{ui,Uj) := (TTi (g) ■Kj){'Rij) G End(Vj g Vj). (4.5) 

For reader’s convenience, we have written it explicitly in the 
appendix. Let P = J2i<a,b<4i~^)’'^a,b g F'ba denote the 
graded permutation operator on the graded vector space. 
We have that R(u, u) = P and the unitarity property now 
reads as R{ui,Uj)PR{uj,Ui)P = I, where / is the identity 
matrix. Moreover, R satisfies the graded Yang-Baxter 
equation in End(I4 g g Vk) 


RijRik^jk — R'ikRikR', 


■jkXlikXlij , 


where Rij = R{ui, uj) g J, Rjk = / g R{uj,Uk) and Rik = 
(J g P)R{ui, Uk){I g P). 

The matrix R{ui,Uj) is the Beisert’s AdS/CFT world- 
sheet R-matrix found in [18]. It is equivalent to Shastry’s 
R-matrix up to a similarity transformation and the follow¬ 
ing identification of the parameters |H] 


xf = I 


a( 0 ) 

b( 0 ) 


, 2 h(e) 


X, = -i 


a{e) 


, 2 h(e) 


(4.6) 


with a^ {0) + b^ (0) = 1 and h a{9) b(0) = 2^ sinh 2h(6*). The 
standard choice is a{9) = cos 9 and b( 6 >) = sin0. 


T,; 






- + h.c.) 


+ s!?(n4n4 -h n4n4)(n4 - n4)(n4 - n^) 

+ - l)(n*T - njt)(n4 - n4) 

+ 2 s|?(niT - 5)(n4 “ l) + “ 5 ) 

+ - l)(n*t - ng)^ (njt - 

X - 5 ) - - 5 ))- (4-2) 

The Boltzmann weights are 


(±) ^ f Vij \ 7i7j 

V1 - xj X- -x+J vfvj ’ 
72 ) - ^7 ) ~ ^7 ^7 ) 73 ) 

^ij / — -t-^/ + — — -t-\ ’ 

(x, -xj)(xl -X, x^ xl) 


(4) 1 

4 = 


„(5) 




-l-l- 
ij J 


3 4 — 


(4.3) 


Lax operator. In this letter we use Latin letters i, j to index 
lattice sites, which are quantum spaces of our integrable 
model. We will use Greek letters A, 7 to index auxiliary 
spaces associated with the model. In our setup, auxiliary 
spaces are isomorphic to the quantum spaces. Thus the 
Lax operator is obtained by associating one leg of the 
oscillator R-matrix with the auxiliary space and mapping 
it to the matrix representation. By this we mean that 

'■= £ End(V)j) g OsCi. (4.7) 

Eor any matrices A and B, introduce a graded tensor 
product by 

{A®B)acM = (-I)“("+4)AbRed. 


This allows us to write the Lax operator in terms of 
the following decomposition (c./. USD 

— K (Zj ^2'!' ^ ^fii H” 2 '^ 4 Q 


where is cl Lax operator for a fermionic lattice model 
with two-dimensional sites 


r . — 

fiia. — 




_p-*0a,V- 


-1 
2 -^ (j-i 


2 )”*” 2 ft 


..ct \ 
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The matrices K^i, P^i and Q^i all have the same form, 
= \{k% + I + i(fc+ - 0 a^ 


For any graded matrix A G Endtdj, introduce a supertrace 
by StrA = X]i<a< 4 (“)’^(^)“a' transfer matrix is de¬ 
fined as the supertrace of the monodromy matrix 


where is the third Pauli matrix. For P^i arrd one 
needs to replace with and respectively. The 
matrix elements in the expressions above contain 




„(2)„(3) 

^ia ’ 


fui — 


ktk- 


and 


<l^i + Pt = (1 - (1 - 

%i - P'ili = (1 - (1 - 

%i+Pli = 4((1 - + (1 - 41 ^)). 

(5)^ 




- Pui = 4(1 - s„ J (1 - 


where 




- 'Suj)('SiL + ■Siu) xt 


fll/ \ IfJ, 

44'(.g^-l) 




xfXi -vl' 


(4.9) 


The inverse of the Lax operator is obtained by a certain in¬ 
version of the parameters x^ of both auxiliary and quantum 
spaces, namely 




(4.10) 


The Lax operator satisfies the fundamental commuta¬ 
tion relation in End(V)j 0 Vx) (S) OsCi 

0 I){I 0 Lx^) = (/ 0 £az)(£m^ ® (4.11) 


which follows from (4.4). The intertwining equation (4.1) 
for bosonic operators in (3.5) now reads 


= [tT;, (3;x ),£/.*]• 


(4.12) 


Notice that the commutator on the left hand side acts on 
the oscillators only, while the commutator on the right 
hand side deals solely with the matrix structure of 

Let Qi = Qiap and be a shorthand notation 

for fermionic operators in ( |3.7[ ). We have the identity 

4“ (Q/i)ll/i2^2; (4.1'1) 


which tells us how to move a local fermionic operator 
through the Lax operator. A similar identity, with and 
Vi replaced with their inverses, holds for Q|. 


Transfer matrix. We introduce the monodromy matrix in 
the usual way 


7^ := G EndV)x 0 Osc 




(4.14) 


It satisfies a commutation relation in End(E^ 0 Ea) <?) 
OsCi...L similar to (4.11) 


<8) /)(/ 0 Ta) = (J 0 Ta)(T^ (g) /)i?^A. (4.15) 


:= Str^T^ G Osci.,. 1 ,. 


(4.16) 


Relation (4.15) implies that [t^,t-a] = 0. Thus generates 
a tower of mutually commutative operators that charac¬ 
terize this integrable model. Note that \ u) is a 

function of the spectral parameters of the auxiliary and all 
the quantum spaces. 

Let us consider a homogeneous lattice. Each quantum 
space is now described by the same spectral parameter 
u = Ui for all 1 < i < L. By expanding the transfer matrix 
in the neighbourhood of = u, we have 


L-l 

— 1 ^ ^ii+i 

i=l 


1 -f (up — u)0~C .. 


(4.17) 


where TC is the Hamiltonian of the general model. We will 
compute the Hamiltonian explicitly in the next subsection. 
We now want to focus on the symmetry properties of 
and 


/i- 


Both and have the usual symmetry properties for 


the bosonic operators in (3.5). They are 

[T^,a] = [7r^(a;.),T^], [a,v] = o. 


(4.18) 


However, for the fermionic operators in (3.7), the symme¬ 
try relations are of non-abelian type that usually appear 
in models with quantum deformed symmetries. For the 
monodromy matrix we find 

v^T^aiQ) = a(Q)T^ - T^7r^(Q^) + (4.19) 


which can be seen from (4.13). This implies that 
Vf,T^a{Q) = cr(Q)T^ - (1 - Str^(7r^(Q^)T^). (4.20) 


The presence of the braiding factors makes the fermionic 
symmetries non manifest. The operator Str^(7r^(Q^)T^) 
is a transfer matrix for a model with twisted boundary 
conditions. In the = 1 limit this term does not con¬ 
tribute, and the fermionic operators relates eigenstates of 
the transfer matrix whose eigenvalues differ by a factor of 
Otherwise, the fermionic operators relate models with 
regular and twisted boundary conditions. We will discuss 
the = 1 limit in the subsections below. 


Hamiltonian. We are interested in the nearest neighbour 
Hamiltonian 94 = ^ IKii+i, which is obtained by taking a 
logarithmic derivative of the transfer matrix. This follows 
directly from the expansion of the transfer matrix and 
definition of the Lax operator. By choosing a convenient 
normalization, the Hamiltonian density is 




Ui—Uj—U 


(4.21) 
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where di = The parameters parametrize the 

Hamiltonian of the general model for a homogeneous lattice. 

A lengthy but direct computation reveals that the 
Hamiltonian density can be elegantly presented if we de¬ 
composed it into four terms. For this purpose we introduce 
functions 


hi 

h± 


— X 

x+ + x~^ 


h 2 — 


1 -h 


(x^) ^ —x^ 
l/±l - I/Tl • 


1 — X x~^ 
x+ + x~ ’ 


(4.22) 


The Hamiltonian density then has the form 

dfy = hi%ij — 2h2'\lij — (/i 2 + 7,7^) + 77 ^■ (4.23) 

Zfli 4/Zi 


The kinetic term 3Cij is 






a^f3=n 

^j /3 = 1 - h±{n,fi + njfi - l)-{h+ + h_){nip-l){njf 3 -l). 
The term is the usual Hubbard model potential 

Vy = (n,t-i)(n 4 -i)-h(njy-i)(nj- 7 -i). (4.24) 


The term Uy has a simple interpretation in terms of the 
spin and charge operators. Set 


§3 = 


§i .2 = + 

~ 7(3^i,3+£i.3)i 


S ch 1 nr 

i,3 — 2 5 

osp _ 1 nr 

^i,3 ~ 


These operators form two copies of the su(2) algebra and, 
for 1 ^ n, 6 ^ 3, satisfy ^ ^ ^ahr.^i.r^ 

where Cabc is the Levi-Civita symbol. Recall that the Hamil¬ 
tonian density of the Heisenberg model, that describes a 
one-dimensional lattice of spin-^ particles with spin-spin 
or charge-charge interactions, is 


K 


ipjch X '' nsp/ch nsp/ch 

■3 2-^ 

l<a<3 


(4.25) 


This allows us to write the term Hij as a difference of the 
charge and spin Hamiltonian densities 

U., = - ‘K^. (4.26) 


The symmetries of the Hamiltonian density follow di¬ 
rectly from the symmetry properties of the oscillator R- 
matrix. For instance, both the charge and spin su(2) sym¬ 
metries are manifest. This is however not true for the 
fermionic symmetries. Let us consider the fermionic inter¬ 
twining symmetry of 3?^^. We have (for clarity, we write 
the Ui and Uj dependence explicitly) 

(Qi(ui) v{uj) IRy (iti, Uj) 

— (^z(^i) 4“ ^(^ 2 ) ^ 7 (^ 7 ))' 

Expanding this identity to first order in the neighbourhood 
oi u = Ui = Uj and multiplying both sides of it by Tij , we 
find 


[:Kij{u),Q3u)+v{u)Qj{u)] 

= (i/(m) - l)(5Qi)(u) - (9i/)(u) Qj{u), 

where d = ^. In other words, we find that even though 
the oscillator R-matrix has fermionic symmetries, they are 
not manifest at the level of the Hamiltonian density. 


Jordan- Wigner transformation. To make the double fer¬ 
mionic nature of the model manifest, it is useful to introduce 
the so-called Jordan-Wigner transformation. Consider two 
families of mutually commuting Pauli matrices, erf and rf 
with a = ±, z and 1 < i < L. The Jordan-Wigner trans¬ 
formation is the map JW : OsCi,,,i, —>■ EndVliL defined by 


C 2 .T ^ > 

C2.t^(ar---a£)(rf...Tti)r-, 


(4.29) 


where V 4 L = {V 2 is a 4L-dimensional vector space 

equipped with a natural left action of the Pauli matrices Ui 
and Ti (on the copy of two-dimensional spaces V 2 and 
Vf in V 4 L, respectively). The map JW for (3^^ is obtained 
by hermitian conjugating (4.29). 

Applying the Jordan-Wigner transformation to the 
Hamiltonian (4.23), we find the Hamiltonian for a coupled 
spin model with kinetic and potential terms given by 

K^i+l = (f7+cr”+iN+_^^^ -h a-at+i^u+ir) + (o’ ^ 

^U+lr = l-7^±('r/+Tf+i)-|(/l+-4/l_)(Tf-l)(Tf+i-l), 


U 


22+1 


_ 1 (I 

— 8 W* A-i -1 + 


T,- cr; 


2 + 1 


- 5 (+cr,, 1 - a, cri i)(t+T, , 1 - Ti T- 


2 + 1 ^ 


'' 2+1 ^2 

+,+i = i(++ + a^,+,), 

where = JW(!JCij) and similarly for the other terms. 


All the terms in the Hamiltonian density are individually 
invariant under the partial particle-hole transformation, 
[P^,%ij] = [Pi,Uij] = [+,Vy] = 0 , except the constant 
term jh/hi, which is anti-invariant. 

Set (fa = 7 r /2 and zh = 4t/. Use the parametrization 
(4.6) together with the identification h(0) = —£ and 9 = g,. 
Then the Hamiltonian density (4.23) coincides with the 
one given by equation 112.2291 in (20|. up to the sign 
of the constant term, or in other words, up to the P^- 
transformation. 


5. A, B and EKS models 

There are four limits of the general model, that corre¬ 
spond to special points on the elliptic torus parametrized 
by Zhukovsky variables x'^ (see explanation in Section 2.2 
in m)- In this section we discuss three limits that corre¬ 
spond to the A and B models of Frolov-Quinn [15] , and the 
Essler-Korepin-Schoutens (EKS) model [16]. We also de¬ 
scribe the Yangian symmetries of the A and B models. The 
fourth case, the Hubbard model limit, will be considered 
separately in Section 6 . 
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The A model. The A model is obtained by setting 


= 0 , 


x+ = — 


(5.1) 


The Hamiltonian parameters (4.22) in this limit become 

hi =i cot 6, h2=icsc9, /i± = l±tan0. (5.2) 

The symmetries of the A model have a non-trivial braiding 
and consequently the fermionic symmetries of the R-matrix 
do not lead to manifest symmetries of the Hamiltonian. 
However, both spin and charge operators are symmetries 
of the Hamiltonian 


[:K,3a]=0, aG{£,T,lK}, a =1,3. 


(5.3) 


As in the Hubbard model, these su(2) symmetries can 
be enhanced to Yangian symmetries. Remarkably, the 
Yangian symmetries of the A model follow from the ones 


of the general model displayed in (3.8) by making use of 


the fermionic operators. One can check that 
af := 1 + 2isin0 + [Ql; 3 ,aa]Qa /3 (5.4) 

a.,j3 

commutes with the Hamiltonian of the A model. Moreover, 
this model also has a secret symmetry given by 

:= (B + 2^(l + cos^ 0) CSC0 ^ (5.5) 

a./3=t4 


EKS model. The Essler-Korepin-Schoutens model is con¬ 
tained in the general model for a specific value of its cou¬ 
pling constants. In particular, this corresponds to the case 
U = h = iJL = Qoi equation (6) in [TB]. This limit is 
obtained by choosing a different parametrization of . 


Instead of (3.2), we need to set 




+ ^/wf-4 


= x{uidsz\h). (5.9) 


Then the above case of the EKS model is obtained by 
setting (pet = ^ and taking the u ^ oo limit. As in the B 
model, we need to normalize the general Hamiltonian by 
an additional power of u~^ to obtain a hnite result. This 
gives 


hi — 0 , h2 — — 


2 > 


hi h^ ^ 2 ; 


2 hi 


= 1, (5.10) 


which means that the kinetic term 3Cij does not vanish, 
rather its terms are now defined by 




b/3 


= ±i(ni^ -h nj0 - 1). 


(5.11) 


For the EKS model, all the fermionic symmetries are man¬ 


ifest since Vi = 1 in this limit (i.e. the last term in (4.20) 


vanishes). This model actually has a manifest u(2|2) sym¬ 
metry. We note that the auxiliary space of the EKS model 
is different from ours (and that of the A and B models). 
One also needs to take a similar limit of the auxiliary 
parameters x'^. 


The B model. The B model corresponds to the case 


= —X = ^ cos 9. 


(5.6) 


The general Hamiltonian diverges in this limit and has to 
be normalized by multiplying by (a:"*" -|-a;“) to yield a finite 
result. By doing so we find 


/ll=2^cos0, /i 2 = sin^0, = — 2sec6*sin"^-. (5.7) 

h 


Note that the term 777 —Wj in (4.23) vanishes in this limit. 


2 hi 


All the symmetries of the general model are present 
in the B model as well. However, the parameters in the 


symmetry generators are different from the ones in (5.6). are 


In particular, we have to choose different values of the 
phase factors and identify x'^ ~ i so that 


= s’’’; 


X = —X = —^ CSC t 


(5.8) 


By evaluating representation parameters in this way, we 


find that Yangian operators in (3.8) are symmetries of the B 


model. Furthermore, we reproduce the fermionic operators 
Qao-j Qaa givcn in Section 2.3 in [15] by setting the phases 
= Itt for Qaa, Qia ^nd = -^TT for Qa^i3, Qa^/ 3 ^ re¬ 
spectively. Finally, the B model exhibits a secret symmetry 


upon setting pa = 0 in (3.9). 


6. The Hubbard model 

In the remaining two sections we consider the Hubbard 
model limit of the general model. For ease of notation, we 
index with cx) all the relevant operators in this limit. 

Lax operator. The Hubbard Lax operator is obtained by 
setting pa = 7r/2 and taking the limit 


:= lim L 

^ 11 -—4 


fll 




( 6 . 1 ) 


of the quantum spectral parameter Ui. In this limit the 
spin chain is homogeneous. We find that lim.i,^_,,oo = 


lim„.^oo q'pi = 0 , thus = 0 and becomes of 

the standard form. The coefficients of the matrix elements 










where we have suppressed the quantum space index i. In the 
standard parametrization (4.6), = e*®, = sec 9 and 


sec 9. The inverse of the Hubbard Lax operator 
is obtained by substituting x'^ —)■ l/xj, or, equivalently, 
by h( 6 i) ^ -h( 6 »). 

The Hubbard Lax operator satisfies the spin and charge 
symmetries (4.12), However, the fermionic symmetry (4.13) 


does not have its counterpart in the Hubbard limit, since 
Vi —>■ 00 in this limit. 
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Transfer matrix. The Hubbard monodromy matrix T“ 
and transfer matrix are obtained by taking the limit 


operators in p(g) specializes to the Uglov-Korepin Yangian 
operators 


(|6.1|) of all quantum spectral parameters ut and setting 

:= Str^T“. ( 6 . 2 ) 


= 7^/2, 


1 “ := lim T, 


lim ( 03 - I 

W—>-oo \ Z 


[33,Qa/3]Ql;3) =-3“- (7-3) 


I0„=7I-/2’ 


Moreover, the specialization 


They satisfy the spin and charge symmetry relations (4.18) 


althought the symmetries (4.19) and (4.20) do not survive. — a{QaBQ^ )+L{u—-)] = -TC (7 4) 

u—fCJO \ 2 ^ j 

^ r. fl—+ I Z 


Hamiltonian. The Hubbard Hamiltonian is obtained by 
taking the u —>■ oo limit of the Hamiltonian density (4.23) 
and setting (j)a = 7r/2, 


:= lim 0f,- 


\ ^Q:=7r/2 ’ 


(6.3) 


In this limit the constants tend to zero, while hi 

tends to unity. Therefore = 1 and we obtain the 

usual Hamiltonian density of the Hubbard model, tK°° = 
3C^ - given by 


= E 

«=n 


,t-V f-i 


), V“=V,,. (6.4) 


Another possibility is to set 4>a = ”77/2, which is equivalent 


to replacing % with —^ in (6.4). This effectively gives a 


Hubbard model with opposite coupling constant, —h. 

7 . A shortcut to the Hubbard Hamiltonian and the 
Uglov—Korepin Yangian 


Yangian operators (3.8) and (3.9) are symmetries of 


the R-matrix, but not of the Hamiltonian (4.23). However 


they become symmetries of the Hubbard Hamiltonian in 
the M —> 00 limit. In particular, spin and charge Yangian 
symmetries specialize to the Uglov-Korepin Yangian ob¬ 
tained in [S], and the secret symmetry specializes to the 
Hubbard Hamiltonian. We will demonstrate this special¬ 
ization for the spin Yangian operators. The specialization 
of the charge Yangian operators is similar. 

Let 0 G {£,T, 5f}. The Uglov-Korepin Yangian opera¬ 
tors in Statement 1) in our notation can be written as 


l<i<L 




2 + 1,3/' 


2 S 

1 < 2 <J<L 


[+3,Z^J], (7.1) 


where the non-local spin symmetry operators are defined 

by 


?U) 

-'i.S 

rU) 

rU) 


- 

— p»(</^4-7-<^T(*+i))f-t r- 

— C C ; I Cj 






(7.2) 


ufU) — I - 


and we have assumed that (j)a = tt/2. A direct computa¬ 
tion then shows that 0§° is a symmetry of the Hubbard 
Hamiltonian on an infinite interval. 

The following linear combination of the spin Yangian 
generators and the quadratic combination of fermionic 


a./3=t,t 

relates the secret symmetry to the Hubbard Hamiltonian. 
Let us now prove these specializations. First, note that 


taking the w —>■ oo limit of (3.7) and keeping the leading 
terms only, we have 

Om = *«•*““ 4. 

aSt = -71 c,i, 

= A-*'*””'.',. 

Upon setting 


a=t.t 






a./3=T,t 


)too _ 

J_p4tv1 


'UT “ 

US® 


II II 

J_e-*04.v. 

Vh 
_ 1 

Vh 

(7.5) 

)too _ 
'4^ “ 

J_ -20.fi 

Vh^ 



'cO-Of-- U 1 

(7.6) 

only. 

we have 


■i+l 


(7.7) 

E 


(7.8) 




Notice that Ca = 2, Cij = —Cji for i ^ j and [2i,3,Cij\ = 
3!+*^ Hence 

H:= y: [+3,w-.+ - y] [a3,Q“^]Qly 

= -l E E +3. 




l<i<L 


In the M —^ oo limit only the j = f -I- 1 terms contribute in 
the sum above. Thus 

f E “32.3+ = - E +3l+y3)> 


, l<i<L 


1<2<L 


which, combined with (|3.8|), gives (|7f3|). 
that 


Next, we demonstrate the specialization (|7.4). Observe 

(7.9) 


E +w+'^) = - E 

E = F E (7.10) 


a./3=t4 
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Therefore 


B:= E 

Oi,0=tA 

= -J E (2-3<5.,>^-^+iC.,. 

l<i<j<L 


and Physical Sciences Research Council (EPSRC) of the 
United Kingdom for the Postdoctoral Fellowship under the 
grant EP/K031805/1 ‘New algebraic structures inspired by 
gauge/gravity dualities ’. 

Appendix A. 


and 

lim E E *(4aC*+l,a +4+l.aC*a) +Ew 

u—^-oo Z ‘ ^ 

l<i<L Q:=t,4' 

which, combined with the identity 


The Cartan matrix and the distinguished Dynkin diagram 
of the centrally extended su(2|2) superalgebra are 


A = 


2 

-1 

0 



0-0-0 (A.i) 

12 3 


lim —;-— 

M->oo ad + be 


- 5)(n4 - 5 ) + 7 


and (3.9) gives (7.4), as required. The key feature of this 


specialization is the braided non-local terms of the Yan- 
gian and secret symmetry operators. In the u —>■ 00 limit 
these non-local terms specialize to the nearest-neighbor 
terms that describe the hopping terms in the Hamiltonian 


or the hopping-like terms in (7.2). The Hamiltonian is 


related to the ‘opposite’ secret symmetry cr('B) because of 
the definition (4.21) and the fact that u) = 


8. Conclusions 

In this letter we have discussed the algebraic structure of 
a general integrable lattice model described by the centrally 
extended su(2|2) superalgebra. This superalgebra is a 
symmetry of Beisert’s AdS/CFT worldsheet R-matrix and 
Shastry’s R-matrix, and gives rise to a one-dimensional 
Hubbard-Shastry type model. Quantum spaces of this 
model are parametrized by quantum spectral parameters 
and are identified with fundamental modules of this model’s 
superalgebra. Despite the fact that the lattice and the R- 
matrix exhibit both bosonic and fermionic symmetries, only 
bosonic symmetries are manifest in the Hamiltonian. The 
fermionic symmetries enter only via certain non-abelian 
identities. 

The general model contains the A and B models intro¬ 
duced in [H], and a special case of the EKS model [T5] as 
certain specializations. Starting with the Yangian symme¬ 
tries of the general model, we found Yangian and secret 
symmetries of the A and B models. 

Finally, we identified the limit in which the general 
model specializes to the one-dimensional Hubbard model. 
We have shown that spin and charge Yangian symmetries 
of the general model in this limit specialize to the Uglov- 
Korepin Yangian operators. Remarkably, we find a purely 
algebraic interpretation of the Hubbard Hamiltonian as a 
specialization of the secret symmetry. 
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The matrix R(ui,Uj) has the form 


R{ui,Uj) = 

= E +(1-4?) E )Eaa(^E,, 

a—1,6—0 a—0, 6—1 

+ (sij — Sjj) ^ {-)°'Eab® Ei)a-\- 44 E Eaa®Ebh 

a+6=l a—b—1 

+ i4f E +4? 4? E ){Eaa®Ebb-Eab®Eba) 

a^b a^b 

d—b—0 0=6=1 

+ (4-+%) E (A.2) 

a=c, b—d 
a^b^c^d 


with the Boltzmann weights given in (4.31. 
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